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Abstract: We present a method of analyzing a set of N time signals /,(?) that consist of local 
measurements of the same physical observable taken at N sequential locations Z, along the length 
of an experimental device. The result is an algorithm for reconstructing an approximation F{z,t) 
of the field f{z,t) in the inaccessible regions between the points of measurement. We also explore 
the conditions needed for this approximation to hold, and test the algorithm under a variety of con- 
ditions. We apply this method to analyze the magnetic field measurements taken on the Compact 
Toroid Injection experiment (CTIX) plasma accelerator; providing a direct means of visualizing 
experimental data, quantifying global properties, and benchmarking simulation. 
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1. Introduction 

Here we present a new method to reconstruct an accurate estimation of the instantaneous spatial 
dependence of an experimentally measurable field, using the signals from a small set of spatially 
distributed time-domain measurements. Such a reconstruction is only possible if the field is being 
transported past the sequence of measurement locations at a speed that is comparable to the rate at 
which field growth or decay is occurring within individual fluid elements. In this paper the term 
"field" will refer to either a variable property of a fluid medium (such as density or pressure), or 
to an independent quantity (such as an electromagnetic field), which may or may not couple to a 
physical fluid. 

We call this method co-moving interpolation because it follows the trajectories of fluid ele- 
ments through the system and performs an interpolation of the field quantity using the measured 
values at the points where a trajectory crosses the probe locations, doing so along the path of each 
trajectory. The algorithm for this method is described in detail in sections || and Q 

The development of this method was motivated by data from the Compact Toroid Injection 
experiment (CTIX), [jl|], [^] a plasma accelerator device that generates high speed magnetized 
plasma rings called compact toroids (v ~ 200 km/s, B ~ 10 kOauss), primarily for the purpose 
of refueling a tokamak fusion reactor. This interpolation method has been particularly helpful in 
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resolving the issue of the magnetic geometry of the compact toroid plasma, which determines how 
the CT will interact with the magnetic field of the reactor as it deposits fuel at its interior. In section 
^ we will present an overview of how this interpolation method is used on CTIX, while in section 
^ we will look at the particular method of velocity estimation that works well for our system. 

This interpolation relies on first making an estimation of the velocity field throughout the 
system, usually based on the apparent time-of-flight kinematics of any traveling pulses. We treat 
the velocity estimation as an independent problem that is not necessary to fully explore in order to 
understand the key features of the co-moving interpolation algorithm. However, in section |6| we 
will examine an analytic model of shock propagation as a way to compare an exact solution against 
three cases of interpolation using different velocity fields. This comparison illustrates an important 
distinction between the fluid velocity field and the optimal reconstruction velocity field. 

Lastly, in section |8| we will examine two different approaches to find upper bounds on the 
intrinsic errors of this method and show how they scale with the system parameters, (e.g., number 
of probes, flow velocity, probe separation). We demonstrate the application of this method to 
facilitate the direct comparison between simulation, analytic models, and experiment. 

2. Sequential measurements of a moving pulse 

The motivation for developing this type of interpolation originates with the task of trying to inter- 
pret experimental time-domain signals. Effective implementations of this algorithm can produce a 
quantitative analysis of the accelerator dynamics that significantly improves upon previous meth- 
ods. (Compare to [|l|], [|], [|]). 

On CTIX our interest is in reconstructing the approximate spatial dependence of the magnetic 
field within the plasma, using data collected from three magnetic field probes that measure the edge 
magnetic field of an accelerated plasma pulse. Sample data is presented below in figure |l| showing 
the time dependence of the magnetic field measured by probes. In Figure ^ we see a schematic of 
the experimental arrangement in which these measurements were taken. 

It is interesting to note that this general concept is also used in video decompression and key 
frame animation when trying to interpolate in time between original frames of a video image 
[^], although the mathematical formalism and resulting algorithm are not applicable to analyzing 
the data from a set of diagnostic probes, and a new formulation is required. 

Any implementation of this method must begin by making an accurate estimate of the fluid 
velocity field for all {z,t) in the system. Velocity estimation is a difficult problem in general, [Q] 
but many experimental systems have some simplifying symmetry or property that allows the flow 
velocity to be inferred from probe measurements. In the case of CTIX, the simplifying property 
used is the approximation of constant acceleration, as discussed in section ^. 

Fluid elements being tracked may not necessarily correspond to any physical fluid with mass 
density. These may correspond to the trajectories inferred from either the group velocity or phase 
velocity of the wave pattern that is being reconstructed. The interpolation algorithm takes this esti- 
mated velocity field as an input and uses it to create fluid element trajectories, on which it computes 
a superposition of the adjacent real data to determine the interpolated value of the observable along 
each trajectory. 
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Figure 1. Example of CTIX magnetic field data, three measurements of the axial component of the edge 
magnetic field B^{t) taken within the CTIX plasma accelerator at positions z = 57, 91, 142 cm. 




Formation Section 

Figure 2. Schematic of the CTIX plasma accelerator showing the coaxial inner and outer electrodes, the 
formation and accelerator circuits, and the location of the magnetic probes. The co-moving interpolation 
method has been helpful in analyzing the magnetic data for this system. 



3. Co-moving interpolation formula 

This key formula defines the interpolated value F{z, t) of the observable f{z,t) at position z between 
the ith and (/+ l)th real probes, using a superposition of neighboring probe signals evaluated at 
retarded and advanced times: 

F{zj)=Wi{z,t)fi{Ti{z,t)) + {l-Wiiz,t))fi+iiTi+i{z,t)), (3.1) 

where w, is a non-negative weighting factor (< 1) that has boundary conditions Wi{Zi,t) = 1, 
Wi{Zj+i,t) =0. The function /,(?) is the time signal of the if /i real probe =f{Zi,t) where Z, is 
the position of the ith probe, / = 1,2, ...A^. The time shift function T,(z,f) records the time that the 
fluid element at {z,t) crosses through the z = Z,- position. For a fluid that is moving in the positive z 
direction, a given fluid element first passes by the lower probe at Z,-, then passes through the value 
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of z under consideration for interpolation, and finally passes the upper probe at Z;+i. For all z in 
the interval [Z,,Z,+i], we have T,(z,f) <t< T,+i(z,?). The time shifts must equal the identity at the 
probe locations, Zi{Zi,t) = t for all t. 

Although in general the weighting function w, can depend on z and t, an efficient time- 
independent method is to interpolate linearly in space. 

Wi{z,t) = {z-Zi)/{Zi+i-Zi). (3.2) 

Alternatively, a linear interpolation in the time dimension could be performed using 

Wi{z,t) = {t-Ti{z,t))/{Xi+i{z,t)-Ti{z,t)). (3.3) 

This is slightly more computationally expensive than the time-independent method due to the look- 
up time of T, (z, compared to the constant Z,-, but it does a better job at matching rates of change of 
observables that do not stay constant along the fluid trajectories. Higher-order polynomial interpo- 
lation along the trajectories is also possible, but may be unnecessarily complicated given the good 
results of the methods that have been implemented so far using simple linear weighting functions. 



4. Trajectory algorithm for T,(z, t) 

For continuous, integrable velocity fields the trajectory of an individual fluid element given by 
z = z^(0 is related to the velocity v(z,f) according to 

^=v(z,(0,0=v,(0. (4.1) 

Here jU is a Lagrangian coordinate that uniquely labels the fluid elements, which can be defined 
as the position of the fluid element at some initial time t^, by z^(fo) = /i. A second order Runge- 
Kutta algorithm ^ is used to compute each trajectory from the given velocity field v(z,f) and 
find the time that it crosses the neighboring probe locations. For each probe position Z, there is a 
distinct sub-domain of grid points = [Z,_i,Z/+i] x ^,Nt\ over which the crossing time Ti{z,t) 
needs to be evaluated. The trajectory waveform is composed of a time-coordinate array {tj} and a 
space-coordinate array {zj}- The magnitude of the timestep \h\ is a small arbitrary constant that is 
fixed before runtime. 



Integration of ( 1 1 ) is accomplished with the following algorithm. 



Outer Loop For / = 1 , 2 , . . .A/^, set Z, equal to the ith real probe location and evaluate the following 
Steps 1-4 for every (z,f) pair in the ith sub-domain H,-. 

Inner Loop: Step 1 Start the trajectory at t^ = t, zo= Z- Set the sign of the timestep h depending 
on the direction toward the probe at Z,-: /j > if z < Z,, /j < if z > Z,. 

Step 2 Find the next trajectory point using the Runge-Kutta midpoint method 

tj+i =tj + h (4.2) 

= zj + hv(^zj + ^v{zj,tj), tj + ^^ (4.3) 



-4- 



step 3 Repeat Step 2 (iterating 7 = 7 + 1) until the interval between zj and zj+i contains the probe 
location Z,. After crossing the probe location, the final point added to the trajectory waveform is 

iZj+i,tj+i). 

Step 4 Find the crossing time by interpolating across the interval [zj,Zj+i], such that T;(z,f) is in 
the interval [tj,tj+i]. For a linear interpolation of the crossing time, set 



Ti{z,t)=tj + h- 



Zi 



(4.4) 



Refer to figure ||for a visual depiction of Steps 1-4. When our implementation of this method is run 
on a 1.8 GHz Pentium 4 computer it yields a runtime of about 4 iJ.s per timestep on an individual 
trajectory calculation. The local truncation error of the trajectory algorithm is 0{h^). 

Once the crossing times have been tabulated in T,(z,f) for / = 1,2, ...N, the weighting functions 



Wi{z,t) can be evaluated, and the co-moving interpolation formula ( pTTI ) can be applied for all 
in the full domain. 



X 




1 \ I I \ \ \ \ \ \ I r 
Z axis ( 1 cm grid) 



Figure 3. The estimated velocity field on the (z,t) domain. The trajectory of a fluid element is shown as the 
grey curve as it passes the probe at position at Z, and then the probe at Z+i . 
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5. Frozen-in approximation 



The concept behind formula (3.1) is the notion that if an ideal fluid has zero substantial derivative 
for a given observable /, 

^ + v^ =0 
dt dz 

then / is purely convected along by the velocity field, or "frozen" into the fluid, and the velocity 
field alone contains all of the dynamical information of the system. The frozen-in approximation 
Fi{z,t) of f{z,t) about the measurement location Z; is simply 

Fiiz,t)=Mzi{z,t)), (5.1) 

where /,(?) = f{Zi,t) and T,(z,?) is the time that fluid trajectories cross Z, as determined from a 



known velocity field. This is exactly the N = I case of the general interpolation formula (3.1 ), with 

Wi{z,t) = 1, for allz,f. 

We find that for ideal fields that are perfectly frozen into the fluid, the reconstructed field can 
be made to agree with the original model field arbitrarily well. However for real fields that grow or 
decay as they convect the N = I interpolation will have an error that can increase without bounds 
as the distance away from the measurement location Z, increases. 

The full interpolation ( ^ ) with multiple probes will improve the reconstruction provided that 
enough probes are used to guarantee that the characteristic time of either growth or decay of the 
field is longer than the transit time between sequential probes. This requires that the growth or 
decay rate be constrained by: 

where (/) is the average value of the field, (v) is the average velocity, and ALp^obe is the distance 
between adjacent probes; for N uniformly spaced probes, ALprobe = L/{N + 1), where L is the total 



length of the system. If (5.2) is not satisfied then more probes are needed to make the co-moving 



interpolation work well on the given system. 

6. Reconstruction of test case: Burgers' shock 

A simple nonlinear one-dimensional model of shock propagation is given by Burgers' equation 
for the characteristic velocity c(z,?) in a viscous fluid, 

dc dc d^c 

^^'Tz = 'd?^ ^'-'^ 

where v is the kinematic viscosity and provides a diffusion mechanism for smoothing out discon- 
tinuities near the shock front. This has a steady-state shock solution of the form 



C{z,t) =C1 + (C2-Cl) 



l+exp(^(z-[/0" ' 



(6.2) 



where the constants c\ < C2 give, respectively, the characteristic speed of wave propagation far 
ahead, and far behind the shock front, while U is the speed of the shock front itself. In the case 
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Figure 4. Burgers' shock wave tests, a) Original shock wave c(z,f) represented in a cyclic grey-scale to 
emphasize gradients (left). Trajectories of fluid velocity are shown (middle). The steady-state shock profile, 
(right) shown at time t ~ 200. b) The co-moving interpolation using v fluid (z,t) as the velocity field results 
in a fan-out backwards in time, shown with N = I (left) and with = 10 (middle). We then compare the 
reconstructed profiles with the original (right), c) Reconstruction with the optimal velocity v(z,f) = U using 
only one probe (A^ = 1). 



considered here U is the average of the two speeds, U = ^{c\ +C2). The shock thickness is deter- 
mined by the diffusive term, with an e-folding thickness Lshoc/i = 4v/(c2 — ci). 

For a density flux given hy pv / = q = Q{p) — vdp/dz, the characteristic velocity is defined as 
c(p) = Q'{p)- In the case of simple model of flux with a parabolic dependence on p, Q{p) = ccp^, 
the characteristic velocity depends linearly on density c{z,t) = 2ap{z,t). The fluid velocity for 
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this system is 



q V dp c V dc 

Vf = - = ap ^ = ^. (6.3) 

^ p pdz 2 cdz 



This model is interesting as it applies to this interpolation, because (|6JJ) equates the substantial 
derivative of the c(z,f) with vd^c/dz^, and so for v 7^ trajectories that follow the motion of wave 
elements will pass through changing values of c{z, t). This means that the field c(z, t) is not "frozen" 
into the velocity field c{z,t). As it flows, it changes at a rate defined by the viscosity and the second 
derivative of c. The same is true for the fluid velocity v j and the density that is transported by it. 

If we use V f{z,t) as the velocity field for the co-moving interpolation for small A/^ ~ 3, the 
reconstructed field compares poorly to the a priori model. The source of the problem is that the 
convergence of the actual fluid trajectories toward the shock front results in a "fan-out" going 
backwards in time, shown in figure ^. This problem is mitigated if more probes are added to lessen 
the span of the interpolation, but the problem is fundamental. 

We can quantitatively compare the shock wave field c{z,t) given by ( |6^ against the results 
of co-moving interpolation for three different cases on the basis of absolute value of relative error 
\R„r\ = \{F{z,t)-c{z,t))/c{z,t)\. 

Case 1. v(z,f) = \f{z,t), N = I: The interpolation yields max \Rerr\ = 2.41 and a mean {\Ren-\) = 
0.143. This large error is understandable given the fact that c{z,t) is not frozen into Vf{z,t). See 
Figure I b) left. 

Case 2. v(z,f) = Vf{z,t), N = 10: With ten probes the result is better with max \Rerr\ = 0.492 and 
a mean {\Rerr\) = 0.032. The error is largest halfway between probes, which results in a oscillatory 
behavior in the reconstructed shock profile, shown in figure ^ b) middle. 

Case 3. v{z,t) = U, N = 1: The interpolation is near-optimal with maximum max|/?e„| = 2.6 x 
10^^^ and a mean of {\Rerr\) = 2.0 x 10^^, which is at the level of the trajectory integration error. 
See Figure ^ c). 

The conclusion to take from this test comparison is that the translational velocity U of the 
steady-state wave pattern results in a better reconstruction than when the actual fluid velocity field is 
used. This distinction is most important when the fluid is not ideal and viscous diffusion is present. 
When we use the optimal reconstruction velocity v{z,t) = U, then the co-moving interpolation is 
optimal even with a single probe measurement (N = 1). In general, if we are trying to reconstruct 
the observable f{z,t) we would get the best results if we could perform the interpolation along 
trajectories that are the level sets of /. In principle this could be done by computing 

, ^ df/dt 
df/dz 

In the real problem we only have the N probe signals to work with, and the apparent velocity that 
can be inferred from them. As we have seen, this apparent velocity of any steady wave pattern 
is very close to the optimal velocity field even if the actual fluid velocity is known by some other 
means. 
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7. Velocity estimation for the CTIX system 

In systems with freely accelerating flows, a velocity estimation method similar to the one used on 
CTIX is likely to have some success. We estimate the flow field on CTIX by tracking one distinct 
feature as it travels down the accelerator, and then apply the same relative kinematics to the rest 
of the fluid elements in the system. The most stable feature on the CTIX waveforms is the back 
edge of CT; see Figure [T[ This is the junction point between the CT and the pushing field, and we 
can reliably define the arrival of this point as being the time when the B, signal of a given probe 
crosses the half-maximum level on the trailing edge. The high accuracy of this method is due to 
the reproducibility of the steep slope that occurs at the CT back edge. The crossing times at which 
the CT back edge passes by the three accelerator probes serve as input data for a simple kinematic 
analysis of the average velocities, and overall acceleration. This yields a constant-acceleration fit 
for the trajectory of the CT back edge. 



For positive values of acceleration there will always be a minimum z-position, Zmin < Zcrit), de- 
fined by 



This is the point on the trajectory where y{t) =0. However, for early times before the fluid element 
has reached Zmin, according to this description, the CT would have a negative velocity which is 
unphysical for our system. Instead we will only apply the uniformly accelerated model to the 
region z > Zmin^ and during times of positive velocity. If we work with a simple model in which 
all the fluid elements have the same acceleration throughout time and space, they will all have the 
same Zmin, but they would pass through it at different times. Based on this, we arrive at a velocity 
field of the form 



Notice that this depends on space but not on time. When Zmin > we need a different way to handle 
the velocity in the region < z < z,„„,. Typically Zmin < 57cm, which is a region where velocity 
measurements are unavailable. One solution is to truncate the domain, and only work where the 
probe data implies a velocity field. This is fine when possible, but certain applications of this 
method require making a velocity estimate over the entire domain, (such as the interpretation of a 
Doppler-shift measurement along an axial chord). 

A minimal compromise is to set the velocity equal to a constant value that matches the acceler- 
ated velocity curve at some point Zfit ~ Zmin- For instance, we have found good results by defining 
the position of the fit Zfit = min(|z„„„, 57cm), then defining the velocity in a piecewise fashion 
according to 



Physically, a moderately large initial velocity is needed to be in agreement with the fact that 
there is a high rate of magnetic flux input from the external circuit, as well as an ongoing formation 
of new plasma by ionization of the steady flux of neutrals from the slowly closing gas valve. 



ZcT{t) =zo + vot + -dJ:'^. 




\{z,t) = v^2a(z-z,„;„) 



(7.1) 




(7.2) 
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The result of the velocity estimation for real CTIX data is shown in the left plot of figure 
and the resulting co-moving interpolation is shown in the right plot of figure ^ the traces represent 
time slices of the field B^{z,t). One useful observation that is apparent with the interpolated signal 
is that the CT is expanding as it travels, an effect that is not obvious looking only at the raw time 
signals. 

8. Error analysis 

We present two complementary methods of finding the order of magnitude of maximum errors for 
this method, in the following two subsections. 

8.1 First order velocity error 

It is reasonable to suppose that for well behaved functions there should be some input velocity field 
that yields the optimal reconstruction of the observable field f{z,t), such that the error \F{z,t) — 
f{z,t) \ is minimized. 

A simple approach is to consider an optimal trajectory Y that is a straight line with velocity 
Vr. A deviation away from T to some other non-optimal trajectory j8, with \ ^ = Vr + 5V, will 
produce a quantifiable error in the interpolation. 

The difference in crossing time for the non-optimal trajectory is 

AL AL AL 5V 

^tn = = —. 

^ Vr v2+Vr5V 

Trajectory j3 will cross the probe positions at T, + 5t, and t,+i — 5t,+i where the displacements 
are given by 

= (1 — H',)Af^ and 5t,+i = w/Af^ (8.1) 



a) CTIX velocity field b) Co-moving interpolation of B^(z,t) 
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axial position z (cm) axial position z (cm) 



Figure 5. Example of co-moving interpolation of magnetic field data, (a) The trajectories of the estimated 
velocity field using time-of-flight kinematics from probe pulses. The image of the interpolated B^{z,t) field 
is overlaid in a grey-scale (white = 0, dark > 0) to illustrate where the CT is in the system. A constant 
velocity field is used before the accelerator current fully turns on at about 6 jJ-s. (b) The reconstructed axial 
dependence of the magnetic field at three moments in time, f = 8.2, 9.6, 1 1 .2jj.s. 
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P: V0=Vr + 5V 



Figure 6. Diagram of optimal (F) and non-optimal (j3) trajectories for derivation of first order velocity error. 



If we evaluate the eiror between the optimum (F), and non-optimal (j3) interpolations using a 



Taylor series for fi{t) and /i+i {t) about the points T,- and t,+i respectively, formula ( pTTI ) becomes 
Fr(z,0 -^>(z,0 = wKl - w,-)A?^ -/i^ifei)) (8.2) 



-A?|((i-w,-)/;"(^0+w,-/;^i(§,-+i))- 



+ 



The second term is a remainder that is second order in A?^ . Equation ( |8.2[ ) holds exactly for some 
constants and such that T,- < <^,- < T,- + 5t,-, and T;+i — 5t,+i < < T;+i. 

We see that for fixed F only the weighting function depends on position w, = w,(z,f), and so 
the interpolation error is proportional to w,(z,f)(l — w,(z,f)), which goes to zero at z = Z,',Z,+i and 
it will have a maximum of 1/4 when w,- = 1 /2 (at or near the midpoint z = {Zi + Z,+i ) /2). 

The assumption of straight line trajectories is not fundamental, this analysis generalizes to 
curved trajectories yielding an identical error bound. The details of this generalization are unnec- 



essary for this paper and contain only notational complications. Equation ( |8.2[ ) holds for any curved 
trajectory, the only change is that the term Atp contains higher order terms in 5V that do not affect 
the order of magnitude of the error for small 5V. The error is bounded by a maximum value, 



max 



Ft -Fa 



O 



ALdY 



O 



LdY 



(8.3) 



where L is the total length of the z domain, A'^ is the number of probes used to take measurements. 



and ^ is the characteristic time defined in (pJ|). The interpolation error is proportional to the local 
velocity error 5V. This result demonstrates the advantage of a high speed flow, and the use of many 
probes. 

The first order term of (p^). 



iFr-Fp)^,)=Wi{l-Wi)Atp{f!{xd-fUri+i)) (8.4) 
can be evaluated for real data, given a velocity field, to yield upper and lower error bounds on a 



reconstruction. To demonstrate the validity of this error bound, the necessity of condition (5.2), and 
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the overall performance of this interpolation we have performed a more difficult, N = 2 interpola- 
tion using only the probe signals B,(57,f) and B^{l42,t) as input data, to attempt to reconstruct the 
actual magnetic measurements from the z = 91cm probe. To optimize the reconstruction we will 
use the estimated velocity field that is implied by all three probe signals (57, 91, 142). 




-1 - 



-2 - 
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time (|ls) 

Figure 7. Test comparison with real data. N = 2 interpolation at z = 91cm ( — Test_Bz_ 91) using only 
B,-(57,f) and B,-(142,f) as input data, compared to the measured probe signal B^{9l,t) ( • ) for the same 
shot. Equation (8.2) is used to estimate the upper and lower bounds of the reconstruction error. 



To apply ( |8.4D to this reconstruction we will need to assume that our velocity field is close to 
optimal. The estimated velocity error 5V has a maximum of 5 cm/ i^s. For our probe locations the 
linear weighting function (3^) yields w,(l — w,) = 0.24. In figure ^ we see the good agreement 



of the N = 2 reconstruction and the real probe signal at z = 91 cm. However, within the precursor 



plasma at times f = 4 to 8 jxs, the growth/decay rate condition ( |5.2| ) is not met by this test system 
with only two probes. During this early time period, the relatively light precursor plasma evolves 
on a faster timescale than the body of the CT. The large precursor of about 4 kGauss at the z = 57 
cm probe decays down to 1.3 kGauss by the z = 91 cm probe and then stays almost constant for the 
remainder of the acceleration, see figure [T[ This demonstrates the critical importance of condition 
( p^ ) in designing a probe array. For the precursor plasma, N = 2 is not enough; the precursor 
decay timescale is 2.7 /is, while the temporal span of the N = 2 interpolation is 4.2 iJ.s. We see that 
A'^ = 3 is good enough, since the span of the interpolation is 2.1 /is. 

Figure also demonstrates that for the main structure of the CT (? > 8 /is), where condition 
( P^ ) does hold, the estimated error ( ^^ ) coiTcctly bounds the actual error in the vicinity of the 
compact torus. Thus, condition (|5.2[) is needed both for an accurate interpolation, and also to make 



sohd error estimates for use in subsequent data analysis of the reconstructed fields. 
8.2 Constraints on the amplitude of undetectable transient pulses 

We will now take a different approach and consider the maximum errors that could exist in the 
space between consecutive probes. This will include errors that are as large as mathematically 
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allowable, with no concern initially for the physical limitations on such errors. The result of this 
analysis is a constraint in the form of an uncertainty relation between the duration and spatial extent 
of any undetected transient pulse-like modulations of the field occurring in the unmeasured region 
between the probes. This will provide a solid upper limit to the total error of the reconstruction. 

Since the probes themselves yield accurate and reproducible measurements of the field quan- 
tities in the immediate vicinity of the probe locations, there is little error due to global fluctuations 
of the real field, since these would be detected simultaneously by multiple probes. Instead, the real 
cause for concern are transient fluctuations, that cause error because they are short enough in dura- 
tion and spatial extent, and happen to occur deep enough into the empty space between probes so 
that they go undetected, and consequently the interpolation method has no ability to include their 
existence in the reconstructed waveform. 
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Figure 8. Transient pulse at various times and its time-envelope as a function of axial position. 



It is informative to consider transient pulses that have a Gaussian form in space and time, 

^{z,t)=Ae e y^^) , (8.5) 



where A is the amplitude of the pulse, T its duration, i the characteristic spatial extent, tp the time 
at which the pulse reaches its peak value, and v is the pulse group velocity, which we will assume is 
approximately the average fluid velocity. For a given probe separation ALprobe it is useful to define 
a dimensionless pulse extent equal to -|- V'^T^)/AL^j.^^^. To measure the error introduced by this 
pulse in our reconstruction, we need the envelope ^emiz) of the pulse as a function of position z, 
which is the maximum value of the pulse at a fixed z, over all values of time. 



,{z) =Aexp - 



(8.6) 



In order for this pulse to be undetected by neighboring probes located at Z, and Z,+i we would 
need the pulse envelope to be less than some noise floor Cnoise 



and 



^env(-Z;+l) < £n 



(8.7) 
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where £„o(.ve is the expectation value of fluctuations within an ensemble of measured signals. This 
will have contributions from electrical noise as well as experimental irreproducibility. It is possible 
that real fluctuations could occur at a larger than normal amplitude and in a transient fashion be- 
tween the probes, and thus go undetected. In the worst case, the peak of the transient pulse would 
occur at exactly the half-way point between consecutive probes, thereby taking advantage of the 
largest possible amplitude allowable by ( |8.7D . In this case, we could make the most conservative es- 
timate for the accuracy of this reconstruction method. The worst case is when \tp = {Zi^\ +Zi) /2, 



and so let AL 



probe 



■■2(Z, 



1+1 



vtn). Then condition for non-detection is 



Aexp 



^^probe 
'4(^2 + v2^2) 



< en 



(8.8) 



The maximum possible undetectable amplitude is constrained by 



A < e„ 



AL 



,exp 



probe 



4(£2 + v2t2) 



(8.9) 



This inequality is plotted in figure ^ where the curve indicates the maximum possible ampli- 
tude for an undetectable pulse of characteristic length i and duration T, given a probe separation of 

probe ■ 



o 2 - 



B 1 



Region applicable 
to CTIX data 



A: 



I 

0.412 



Pulse Extent ( h v\^)/AL^ 



probe 



Figure 9. The maximum possible amplitude of undetected transient pulses, plotted as a function of pulse 
extent. The maximum amplitude is determined by the inequality (8.9). 



We can use this result to examine the validity of co-moving interpolation being applied to 
the CTIX data. For an average fluid velocity of v = 20 cm/ jxs, in the region between Z,- = 91 cm 
and Z(+i = 142 cm, if there exists an undetected transient pulse that spans the space between the 
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probes, and it exists for at least the transit time between probes 2i = 5\ cm ^ £ = 20.5 cm and 
2t = (5 1 /2Q)}jLs =^ T = 1 .215}JLs (see figure then the pulse extent is 0.412 and the amplitude of 
the pulse can be no larger than 

^max — 1.84 ■ Enoise- 

For pulses significantly shorter in duration and extent, the upper limit on amplitude becomes 



much larger than this example. In fact, formula (8.9) goes to oo in the limit of zero pulse length and 



duration. However, since such narrow fluctuations, if they exist, should occasionally occur directly 
at a probe location, their absence at the probes implies their absence between the probes as well. 

If there is some physical reason for some finite upper bound on A, such as conservation of 
energy perhaps, then we can restate equation (E^ in the form of an uncertainty relation 



A\n{A/Enoise) 

For some fixed A > Enoise^ there is a trade-off between the spatial extent and the duration 
of undetectable pulses, and that this relationship depends on the fluid velocity and the distance 



between consecutive probes. In agreement with the result from section |8.1| , inequality ( |8.10| ) shows 
that the best reconstruction occurs when consecutive probes are closely spaced, and there is a high 
flow velocity. 

9. Conclusion 

Co-moving interpolation is a method that is applicable to many experimental and industrial prob- 
lems where the properties of a fluid flow are being measured at several discrete locations, and the 
global dynamics of the system need to be determined. The method is straightforwardly applied, 
and makes efficient use of modem computer power. High accuracy can be achieved with a minimal 



number of probe points if conditions such as ( p.2[ ), ( |8.3D and ( |8.10| ) are employed in the design of 
diagnostic systems. 

We have presented the details of a efficient algorithm for this method, which we have imple- 
mented on a desktop computer. We have run our method on a sequence of test problems in which 
an analytic model of shock propagation is compared to reconstructions by co-moving interpolation. 
These tests demonstrated the need for an optimal velocity field, its distinction from the fluid veloc- 
ity field when diffusion is present, as well as the high degree of accuracy that is possible when the 
optimal velocity field is used. 

We have applied our method on the CTIX system, and presented the resulting co-moving in- 
terpolation of B^{z,t). In future articles we will apply this method extensively to the CTIX plasma 
system to investigate MHD effects during the acceleration process. [^] In addition, the results of 
this method are allowing a new kind of comparison between experimental data and ongoing MHD 
simulations of the CTIX system. Ordinarily, multi-dimensional computational studies of physical 
systems produce results that can be directly visualized, or to make comparisons to experiment the 
simulations can be sub-sampled to produce virtual diagnostic signals, which are directly compared 
to raw diagnostic signals from the corresponding experiment. An alternative method of comparison 
is to start with a set of raw experimental data from multiple diagnostics, and then apply a method 
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of analysis to synthesize the approximate global behavior of the system, and then make the com- 
parison to simulation. Both routes have advantages, however, at present there are few tools to solve 
the inverse problem on the experimental data to make use of the second route of comparison. The 
combination of animated visualization of experimental probe data, and quantitative analysis on the 
reconstructed fields, is providing new insight into the nature of the magnetized plasma in CTIX. 

Lastly, we considered two methods for making error bounds on this interpolation. In the first 
analysis we found good agreement of the first order velocity error ( |8.4[ ) with the actual error of 
the N = 2 reconstruction of B^(91,f) signal from CTIX. The second analysis showed that for the 
conditions present on CTIX, undetectable transient pulses of significant size and duration would 
be limited to a level that was within a factor of 2 of the noise floor of the probes. Both methods 
derived similar relationships between flow velocity, probe separation and the resulting error of the 
interpolation. 
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